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The problem of premixed flame propagation in wide horizontal tubes is revisited. Employing 
the on-shell description of flames with arbitrary gas expansion, a nonlinear second-order differential 
equation for the front position of steady flame is derived. Solutions to this equation, obtained 
numerically, reveal two distinct physical regimes of laminar flame propagation controlled by the 
strong baroclinic effect. They differ by the front shape and flame speed, the ratio of the total 
. consumption rates in the two regimes being 1.4 to 1.8, depending on the value of the gas expansion 

7-H ' coefficient. Comparison with the existing experimental data on methane-air flames is made, and 

explanation of the main trends in the observed flame behavior is given. It is shown, in particular, that 
the faster (slower) regime of combustion is realized in mixtures close to (far from) the stoichiometric 
^ 1 composition, with pronounced changeover in between. 
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All laboratory flames are subject to the terrestrial gravitational field whose influence on the process of de- 
flagration is often prominent. In fact, for many premixed flames evolving spontaneously this influence is dom- 
inant as regards formation of the large-scale flame structure. This is because deflagrations involve flows of 
markedly different densities, characterized by essentially subsonic fluid velocities. A simple dimensional analysis 
suggests that under these conditions, the relative gravity impact can be estimated as (acceleration of gravity) x 
(characteristic length) / (characteristic velocity) 2 . It follows that propagation of all hydrocarbon-air flames in tubes 
wider than 1 cm is strongly affected by the gravitational field, the impact being especially pronounced for flames with 
low local consumption rates. Observations fully confirm this conclusion. Under the action of gravity, the front area 
O I considerably increases, and proportionally increases the flame speed. In the regime of uniform movement, i.e., when 
l— "~ '■ the flame speed remains grossly constant, it may be tens times the normal flame speed. For instance, the normal 
flame speed in the mixture of 6% methane with air under standard conditions is about 0.1 m/s, while its propagation 
j> i speed in tubes of lm diameter is 1.5 m/s. 
I/"") ' A vulgar description of the gravitational influence on flame dynamics is that it is a kind of buoyancy effect, and in 
the case of flame propagation in horizontal tubes goes as follows [1] . The light products of combustion moving upward 
tend to occupy space above the heavy fresh gas, which causes the flame front to tilt and spread along the tube. This 
view is actually inadequate, if only because it implies that the process of flame spreading ought to last indefinitely. 

■ Indeed, if buoyancy were the driving force, any tilted flame configuration would be unstable against further front 
flattening, which would result in unbounded growth of the flame speed. As a matter of fact, no description of flame 
behavior in terms of forces can be adequate, because this behavior is essentially determined by the properties of 
the mass flow through the flame front, and it is flow velocities that are the right terms to describe these properties. 
The true mechanism governing flame evolution in the strong gravitational field is the baroclinic effect - the gravity- 
induced generation of vorticity inside the flame front. It is thus vortical motion of the burnt gases that is responsible 
for the formation of titled flame configuration. Flame spreading stabilizes when vorticity generated by the flame front 

■ curvature becomes comparable to the baroclinic vorticity (Cf. Sec. Ill B"2|) . 

In view of practical and theoretical importance of horizontal flame propagation, it was thoroughly investigated 
already at the dawn of modern combustion science. The very existence of the regime of uniform movement was 
established by Mallard and Le Chatelier 0], after which it was scrutinized by Mason and Wheeler Mason 
Coward and Greenwald [f| , and in the famous works by Coward and Hartwell @, 0] ■ The latter summarized and 
extended the results of the preceding works to infer normal flame speeds in mixtures of various composition, and to 
establish relationship between the propagation speed and tube diameter, which was shown to be more complicated 
than the simple power law. 

Despite extensive experimental data acquired on the horizontal flame propagation, there has been lack of adequate 
theoretical description of this process. The reason for this is obvious: The large increase of flame velocity caused by 
the gravitational field implies that the process is strongly nonlinear, and hence is out of reach of the classic approach 
based on explicit solving of the bulk flow equations. The necessity to solve these equations explicitly is indeed the 
stumbling block of this approach already in the case of freely propagating flames, which restricts its applicability to 
flames with weak gas expansion. Namely, the fresh-to-burnt gas-density ratio, 9, must be close to unity in order to 
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damp the nonlinear effects [8|. Under this condition, the horizontal flame propagation was considered in [9(- The 
relative flame velocity increase within the small (8 — 1 )-cxpansion is small by definition, as is vorticity produced by the 
flame. At the same time, laboratory flames have 9 — 5 — f 0, and produce strongly nonlinear rotational flows of burnt 
matter. Therefore, application of the weak gas expansion approximation to flames in horizontal tubes can be only of 
academic interest. Furthermore, the above remark on the role of the baroclinic effect shows that the potential flow 
models such as Frankel's [l(|, or the bubble model [TTJ] are also irrelevant. The latter, for instance, replaces the flame 
by an open inert bubble moving in a horizontal channel, and gives for the steady flame speed the value OA^^/gb (g is 
the gravity acceleration, b the channel width), which badly disagrees with observations. As was already mentioned, 
relation between the flame speed and tube diameter (channel width) is not that simple. 

A way round the problem of solving the bulk equations was found in [l2|, [l3j where exact equations for the fresh 
gas velocity at the flame front (the on-shell velocity, for brevity) were obtained for two-dimensional (2D) steady flame 
with arbitrary gas expansion. Specifically, it was shown that the velocity field of burnt gases can be decomposed into 
rotational and potential components in such a way that the on-shell value of the former can be explicitly expressed 
in terms of the on-shell fresh gas velocity, while retaining incompressibility and boundedness of the latter. These 
properties of the unknown potential component were then used to exclude it from the set of jump conditions at the 
front by applying an integral identity (the dispersion relation) following from the Cauchy theorem. The result is an 
exact complex integro-differential equation (the master equation) for the on-shell fresh gas velocity and the flame 
front position, which together with the evolution equation constitutes a closed system. This on-shell formulation 
was subsequently generalized to unsteady flames [l4l Il5j , and applied to various essentially nonperturbative problems 
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The purpose of this paper is to give rigorous theoretical description of horizontal flame propagation. Using the 
on-shell formulation, it will be shown that the problem of uniform flame movement admits complete solution in the 
2D case. Namely, the set of the on-shell equations will be reduced in Sec. Ill B 21 to a single second-order differential 
equation for the front position. This reduction relies on the large front-slope asymptotic expansion of the master 
equation, constructed in [l6| . Experiments show that in sufficiently wide tubes, flames are subject to a kind of 
turbulent motion observed as a swirling of the flame front. Occurrence of a small-scale structure generally makes 
the large-slope approximation inapplicable. An auxiliary but important result proved in Sec. IIIII is that regarding 
the large-scale flame evolution, the leading-order effect of this structure is to renormalize the normal flame speed; 
the proof follows closely This theorem makes it possible to apply to wrinkled flames the results obtained for 

smooth patterns. Numerical solutions of the second-order equation are found in Sees. Ill B~3l IIV Cl where it is shown 
that for each set of physical parameters, there are two distinct regimes of laminar flame propagation, characterized 
by significantly different flame speeds. A detailed comparison of the obtained results with the experimental data on 
methane-air flames is carried out in Sec. IIVDI The arguments justifying application of the 2D picture to real flames 
are summarized in Sec. IIV Al The effect of heat losses to the tube walls is discussed in Sec. IIV Bl where it is shown 
that because of the large flame spread along the tube, these losses affect horizontal flame propagation in tubes of 
any diameter. After assessment of the calculational accuracy, given in Sec. IIV C 11 theoretical results are used in 
Sec. IIVDI to explain the measured flame speed versus tube diameter in mixtures of various methane concentration. 
It is demonstrated, in particular, that mixtures close to (far from) the stoichiometric composition favor the faster 
(slower) regime of flame propagation. Conclusions are drawn in Sec. [V] 

II. ON-SHELL DESCRIPTION OF HORIZONTALLY PROPAGATING FLAMES 

A. Master equation 

Consider a 2D flame propagating in an initially quiescent gaseous mixture filling horizontal channel of width b with 
ideal walls. Relevance of this picture to the three-dimensional laboratory flames will be discussed in Sec. IIV Al We 
choose Cartesian coordinates (xi, X2) = (x, y) so that the y-axis is along the upper channel wall, y — —00 being in the 
fresh gas. Then the flame front position at time t can be described by an equation y = f(x,t). Let (^1,^2) = (w,u) 
denote Cartesian components of the gas velocity vector v. Even though the flame speed in horizontal tubes often 
considerably exceeds the normal speed, gas flows induced by the propagating flame remain essentially subsonic. This 
allows the fluid be considered as incompressible, and its density, p, be taken constant on both sides of the flame front. 
The fresh gas will be assumed throughout to be of unit density, so that the burnt gas density is 1 /0. Also, it is often 
convenient to use natural units in which b = Uf = 1, where Uf is the planar flame speed relative to the fresh gas (see 
also Sec. [mC]). 

Below, we will be primarily concerned with the steady regime of flame propagation, which is closely related to the 
notion of uniform flame movement. It will be understood that the steady system is considered in the reference frame 
of the flame front. All functions involved are then time- independent, and the argument t can be omitted. In the 
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absence of external gravitational field, the on-shell fresh gas velocity satisfies the following complex integro-differential 
equation [l2l. [l3j 

2^ + (i + ^){h'-^^ ± }=0, (1) 

where ui = u + iw is the complex velocity, [cj] = lu + — cj_ its jump across the front; the subscript — (+) denotes 
restriction to the front of a function defined upstream (downstream) of the front, e.g., u>_(x) = w(x,f(x) — 0); the 
prime denotes ^-differentiation, N = yl + I' 2 j yn = (v-n) is the normal gas velocity (the unit vector n normal to the 
front points to the burnt gas), a = du/dx — dw/dy is the vorticity; finally, the operator 'K is defined on 26-periodic 
functions by 

+b 

(%a){x)= 1 + l ^ x) fdnaWcot^-x + Hm-m})}, (2) 

-b 

the slash denoting the principal value of integral. 

Derivation of Eq. (JTJ is based on a periodic continuation of the flame pattern along the x-axis, which proceeds in 
two steps: first, using the boundary conditions w(0,y) = f'(0) = (which express impermeability and ideality of the 
upper wall), the given flame pattern is carried from the physical domain x £ [0,6] over to x G [—6,0] by reflection 
with respect to the y-axis. The resulting symmetric pattern satisfying 

f(x)=f(-x), w{x,y) = -w(-x,y), u(x,y) = u(-x,y) (3) 

is then periodically continued along the whole x-axis (using the same boundary conditions at the other wall). This 
trick is very convenient as it allows one to avoid explicit account of the boundary conditions, but it is physically 
inadequate in the presence of gravitational field perpendicular to the walls. This is because reflection performed at 
the first step leads to a discontinuity in the gravitational field - its strength, g, being directed downward in the 
physical domain, turns out to be directed upward in the image domain, Fig. [Ha) ■ I n itself, this discontinuity is not a 
problem, because whether or not given configuration is physically relevant depends on its principal feasibility, and in 
the present case it is not difficult to find a setting realizing physical conditions of the symmetric configuration. Namely, 
occurrence of a discontinuity in g suggests that the terrestrial gravitational field is to be replaced in the symmetric 
case by the field of a thin dense layer located near the plane x = 0. Taking its surface density equal to gjl-KG (g = \g\, 
G is the Newton gravitational constant) then yields gravitational field of the same strength as terrestrial, but opposite 
to it in the physical domain, so that the proper symmetric pattern is as shown in Fig. [lib). In this setting, the 
true flame configuration is obtained by a vertical translation from the image domain, (—6, 0) 3 x — > x + b, or by 
flipping the physical domain, (0, b) 3 x — > b — x. Moreover, one can avoid this technical complication by introducing 
an electrically charged layer instead of the gravitating one, and assuming that the fresh and burnt gases are uniformly 
charged like the layer, with vanishingly small (to neglect electric interaction of different gas elements) charge densities 
standing in the same ratio as the mass densities, viz.. 9. But in any case, the symmetric configuration cannot be 
periodically continued in the x-direction, because of interference of the fields produced by the layers: Since the field 
of each additional layer is homogeneous, it cannot be neglected in the physical domain, which makes the continuation 
meaningless. 

Fortunately, the on-shell description can be reformulated in a way that permits application to symmetric flame 
configurations independently of whether they admit periodic continuation. This reformulation was developed in (l8j 
for flames propagating in symmetric channels with arbitrarily curved walls. It uses the method of Green functions to 
explicitly take into account boundary conditions at the walls, which leads to the following modification of the master 
equation 

^ + ('+*){m'-^ + ^/>^} = °. 

where M = Nv™a + . This equation differs from Eq. (JXJ) by the last term in the braces, which is proportional to the 
spatial average of the second term across the channel. As was noticed in the structure of the new term suggests 
that it describes an effective gravitational field along the channel, which appears in the system as a result of the curved 
flame evolution. For steady flame propagation in straight channels, such effective gravity presumably should vanish, 
for the opposite would mean that the flame is self-accelerating. But it turns out difficult to prove this directly, i.e., 
to derive this vanishing from Eq. ([TJ). As to the real part of f^di]Muj + /v 2 _, it does vanish as the average of an odd 



4 



periodic function. However, its imaginary part is even, and there is no apparent reason why it should be zero. Using 
the small (6 — l)-cxpansion, it can be shown that this quantity is indeed zero at least to the fourth post-Sivashinsky 
approximation (which corresponds to retaining terms of the seventh order in (0—1) in the master equation) (2Q| |. but 
the general proof is still lacking. We return to this issue in Sees. IIV C[ IVl 
Equation ((4]) will be used below in the study of uniform flame movement. 

B. Large-slope expansion 

It was already mentioned that the speed of uniform movement in wide horizontal tubes largely exceeds the normal 
flame speed. Experiments indicate that flames propagating in tubes of diameter d < 10 cm have fairly smooth 
elongated fronts, in which case flame speed is directly proportional to the front area. In wider tubes, flames develop 
complex unsteady small-scale structures observed as a widening of the front. The envelope surface retains highly 
elongated shape, but now part of the velocity increase is due to the small-scale wrinkling, which according to Q adds 
40 — 50% to the speed of uniform movement in comparison with the smooth patterns. 

It follows that in the corresponding 2D case of flame propagation in channels of moderate 1 width (b <10cm), the 
front slope is large, |/'| ^> 1, everywhere across the channel except narrow regions near the walls, where the front 
slope drops down to zero. The width of these regions is of order of the cutoff wavelength, A c , related to the finite 
front-thickness effects which determine flame structure near the front end-points. A c is generally larger for slower 
flames, but in no case does it exceed a few millimeters, whereas b 5^2.5 cm in the experiments considered below. Thus, 
flame propagation in channels of moderate width can be adequately studied within the large-slope approximation. In 
wider channels (b >10cm), condition \ f'\ 3> 1 is violated by the small-scale wrinkling, remaining valid only on average. 
As will be proved in Sec. IIII1 the results obtained for smooth fronts hold true for wrinkled flames described in terms 
of averaged flow variables, which makes the large-slope approximation universally applicable to flame propagation in 
channels of any width. 

A fact of fundamental importance for implementation of the large-slope approximation is the locality of steady 
flame patterns at large flow velocities, first discovered in the study of anchored flames This locality means that 
in regions where |/'| ^> 1, front curvature at any point is uniquely determined by the front position and its slope at the 
same point. A similar statement is true of the on-shell gas velocity. In other words, these functions satisfy ordinary 
differential equations. 2 The underlying mathematical reason is that the operator !K , which mainly determines the 
nonlocal structure of the master equation, becomes localizable when the front slope is large. Namely, its large-slope 
expansion reads, in the natural units (we consider henceforth the symmetric pattern shown in Fig. [lja)) 

l 

(5Ca) (*) = {f{x) -i) Jd v a{T]) + 2 a( ~ ??) X (r, - \x\) + ia(-x)(2\x\ - 1) + O (I) , (5) 
where x( x ) is the sign function, 



+1, x > 0, 
-1, x < 0. 



It is because of the step-like asymptotic form of the integral kernel of 3i that the master equation reduces eventually 
to an ordinary differential equation. The other possible source of non- locality is the last term in Eq. Q, which 
depends on the values of f'(x),u-(x),W-(x) at every point of the front, including the regions where the large-slope 
approximation is not valid. A remarkable fact to be proved below is that this term does not prevent localization, 
falling off eventually from the master equation. 

The following version of Eq. ([5]), expressing the action of "K on a derivative, is also useful: 

(ft a!) (x) = (f(x) - i) {a(-\x\) - a(\x\)} + ia'(-x)(2\x\ - 1) + O (~Y (6) 



1 Prom the standpoint of dimensionless criterion gb/Up all channels with b > U^/g ought to be qualified as wide. In the present 
investigation, however, it is more convenient to associate this term with the onset of instability of smooth flame configuration. 

2 In this respect, the term "locality" should not be understood as independence of distant flame elements of each other. On the contrary, 
the asserted locality implies that the whole flame pattern can be restored (by solving an ordinary differential equation) from the 
knowledge of / and /' at any point. In this way, a nonlocal impact of flame anchors was described in [id . 
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where the prime denotes derivative of the function with respect to its argument, a'(y) = da(y)/dy. An important 
result proved in [l6| is that this formula holds true for functions allowed to have a discontinuity at the origin. 



1. Evolution equation and dispersion relation. Boundary conditions 

The local rate of consumption of fresh gas is determined by the so-called evolution equation which expresses the 
normal fresh gas velocity relative to the front as a local functional of the front slope and the on-shell gas velocity. In 
view of this locality, the evolution equation for zero-thickness flames is that for planar flames, viz., — f/N = Uf 
(dot denotes time differentiation), whereas account of the transport processes inside the front adds to the right-hand 
side terms proportional to the front thickness. The relative value of these terms is of order A c /A, where A is the 
wavelength of the front wrinkling. Therefore, they do not affect flame propagation in channels of moderate width. On 
the other hand, it was already mentioned that the results obtained for smooth patterns apply also to wrinkled flames 
upon averaging over the small-scale structure. Hence, it will be sufficient to use the master equation along with the 
evolution equation for zero-thickness steady flames, which in terms of the Cartesian velocity components reads, in the 
natural units, 

u- - f'w- = N. (7) 

In the leading order of the large-slope expansion, it can be written as 

u- = /'(«;_ + 1). (8) 

This is for /' > 0, i.e., x > in the configuration of Fig. HJa); hereon, whenever a local relation appears, it applies to 
the physical domain. 

Let U denote the steady flame speed with respect to the channel. In the rest frame of the flame front, U is the 
velocity of incoming fresh gas. Hence, the longitudinal component of the on-shcll gas velocity, zt_, is 0{U). Next, the 
total mass conservation implies that 

l l 
U = f d v ^lTf^v n _ » I dr)f = /(l) - /(0) . (9) 
o o 

Therefore, of the same order is the (smooth) front slope, /' = 0(U). Equation (|8|) then implies that the other 
(transversal) velocity component w- — 0(1). An equation relating the three functions can be obtained from the 
following dispersion relation 13] 



(w_)' = 0, (10) 



which is a consequence of potentiality, incompressibility and boundedness of the flow upstream. This equation is 
actually contained in Eq. ((4]), which can be seen by applying the operator (1 — iJC) to its left-hand side and using 

i 2 i 1 

the identity JC = — 1 [13]. But within the large-slope expansion employing expression <j5j) or ((6]), Eq. (fT0|) must be 
used separately. The point is that the expansions of M-operator, given in the preceding section, allow only extraction 

of the leading terms in Eqs. ((4]), (|10[) . This is directly related to the fact that the exact identity "K = — 1 cannot be 
verified using these expansions, for the composition of the leading (integral) term in Eq. ([5]) with the undetermined 
remainder is 0(1). 

To extract the leading contribution from Eq. (p~0|) . it is to be noted that the action of Ji on a flow variable essentially 
depends on its parity properties. Namely, Eq. © shows that !Ka' = 0(1) if the function a(x) is even (assuming that 
a = 0(1)), and Dia' = 0(U), if it is odd. Recalling Eq. © and the above estimates for u_,w_, one readily sees that 
the leading contribution to the left-hand side of Eq. (jlOl) is real. We thus obtain the following relation (x > 0) 

f'(x)w-(x) = (l-x)u'_{x). (11) 

Its validity is conditioned by that of Eq. ©, which requires the function f(x) be monotonic in the physical domain, 
and have large slope. As we know, the latter does not hold near channel walls. This is not actually a problem near 
the wall x = 1, for /' = is not a fundamental requirement in the case of zero-thickness flames. It appears only 
upon account of the finite front-thickness effects, which raises differential order of the equation for the front position, 
calling thereby for an additional boundary condition. Note that Eq. (|11[) is satisfied at the point x = 1 even though 
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/' (1) 7^ 0; by virtue of the boundary condition W-(l) — 0. It is therefore legitimate to discard the narrow region 
where /' drops down to zero, and to use Eq. ([lip also near the wall x = 1, which will be assumed henceforth. Things 
are different, however, near the wall x — 0. Although conditions /'(0) = and /'(l) = have the same status from 
the point of view of their genesis just described, differentiability of the flow variables everywhere within the extended 
channel x € (—1, +1) is essentially used in the derivation of Eq. (|10p and of the master equation. In fact, it is existence 
of the short-wavelength cutoff, A c , that guarantees analyticity of the complex velocity, used to derive the dispersion 
relation. This analyticity would be lost if the front had an edge, which necessarily occurs at x = if /'(0) ^ 0. 

Thus, self-consistency of the approach requires vanishing of odd functions at the origin: /'(0) = w_(0) = u'_(Q) = 0. 
But the narrow region near x = 0, where these functions rapidly turn into zero, is still of no interest. To avoid the 
necessity of explicit consideration of this region, we shall proceed as follows. We note that in view of locality of the 
steady flame structure, announced in the preceding section, initial conditions for the ordinary differential equations 
satisfied by f'{x), U—(x), W-(x) can be specified at any internal point of the physical domain where the front slope is 
large, say, x = 0.5. The values of the slope and the on-shell velocity at this point determine the whole flame pattern. 3 
Let us now think of the functions f'(x), u~(x), W-(x) as satisfying Eqs. ((HJ, (|TT|) for all x € (0,1), including the 
regions near the walls. In other words, we use equations obtained within the large-slope expansion to continue their 
solution to the whole physical domain. It will be seen in the next section that the continuation is unique by virtue 
of the uniqueness theorem for ordinary differential equation resolved with respect to the higher derivative. Thus, 
the original configuration and the continued one coincide everywhere except in a vicinity of a; = 0, where the true 
functions f'(x) and W- (x) rapidly turn into zero, but the continued ones do not. To avoid proliferation of symbols, the 
current notation will be reserved for the functions continued in the way just described (in what follows, the original 
functions appear again only once - in formulas p^|) to (|21[) . so no confusion should arise). 

As said, any point with a large front slope can be used for specifying initial conditions, but the end-points x = 
and x = 1 are, of course, most suitable, and they are now admissible for this purpose, as the continued flame pattern 
is local in the whole physical domain. Let wq and /g denote the limiting values for x — > + of the functions u>_ (x) 
and f'(x), respectively: 4 

w o = w_(0+), / = /'(0+). 

These parameters, introduced for later convenience, must be expressed eventually in terms of the physical parameters 
of the theory - 9,U,g. But one relation can be inferred already now from the evolution equation. To this end, we 
have to identify boundary conditions for the functions f(x),U-(x). First, it is a mere convention to set 

/(0)=0, (12) 

but in conjunction with Eq. ^ it implies that 

/(I) = U. (13) 

The last condition follows from potentiality of the upstream flow: By virtue of the Thomson theorem, one has a = 
there, which implies that 

du/dx = dw/dy = 0(1). 
Hence, to the leading order u(x, y) = U for y ^ 0, in particular, u(0, 0) = U, or 

U-(0) = U. (14) 
Now, setting x — in Eq. ([5]), the latter condition gives 

We mention finally that there is a simple algebraic relation between the functions f(x) and U-(x) [ltSf . Namely, 
eliminating u>_ from Eqs. ©, (fTTj) gives u_ = (1 — x)u'_ + /', or [(1 — x")uJ\' + f = . Integration of this equation 
under condition (IT3l) yields 

f = U-(l- x)u_ . (16) 



3 It should be stressed that in general, these values are functionals of the entire flame structure, and in particular, of the flame properties 
near the walls. But it will be shown later on that within the approximations used, they are actually independent of the near-the-wall 
flame behavior (see Sec. Ill B~2t . 

4 The same notation wq appears also in the study of anchored flames |lal , but with a different meaning: wo denotes presently the boundary 
value of the continued velocity, whereas in the case of a flame anchored by a rod of vanishingly small radius it represents the true gas 
velocity induced by the rod. 
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2. Master equation in the leading order. Reduction to ordinary differential equation 



Let us turn to the master equation. The velocity jumps across the front and the on-shell vorticity, entering this 
equation, read in the natural units 
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(17) 
(18) 



The last term in Eq. (fT8|) describes the baroclinic effect [21] . As in the preceding section, the order-of- magnitude 
estimates for the flow variables together with their parity properties imply that the leading contribution to the 
left-hand side of Eq. (j4j) is real. Extracting it with the help of Eq. (|5|), we find 



l r 

2u'_(x)(l + a\x\) ^f- - *yf(x) + f{x) [do ^ ± + a( 



'L 

N 



(v)x(v - M) = o, 



where a = 9 — 1, and 7 denotes the average value of Mw+/v\ across the channel, 

1 

M(o)w+{o) 
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drj- 



v+(o) 



(19) 



(20) 



Just like Eq. (TTT]) . Eq. (TT9")) is valid only for {x : \f'(x)\ >• 1}. By this reason, the large-slope expansion is not 
applicable to the integrand, as it involves the region near 77 = 0. First of all, let us evaluate the integral term with 
(/' /N)'. In the most part of the channel the slope is large, so that f'/N ~ 1 in the physical domain, and f'/N ~ — 1 
in the image domain. Therefore, the function (/' /N)'(x) is localized near x = (as discussed in Sec. Ill B II the slope 
can be consistently assumed large at the wall x — 1, despite the actual flattening caused by the finite front-thickness 
effects). Using this, and taking into account that (f'/N)'{x) is an even function, one can write 



N) 



1 f 

2 N 



= -1. 



Equation (flU)) thus takes the form 



2u'_(x)(l + a\x\) - 



2agx 
u_ (x) 



Mw 4 



af'(x) - 7 f(x) + f(x) / d»7— X (V ~ = . 



(21) 



Dividing it by /', and differentiating with respect to x leads to an ordinary differential equation (x > 0) 



d_ 

dx 



— (1 + ax) 



agx 



f'u- 



a I %- ) (w- — a) = . 

1 u_ u ' 



(22) 



This proves the asserted locality of the steady flame structure, since together with Eqs. (jHJ and (fTTj) . Eq. (|22|) 
constitutes a closed system of ordinary differential equations for the functions f(x), U—(x), W-(x). 

As derived, Eqs. (| 19[) — (f22]) apply for {x : \f'(x)\ 3> 1}, and involve functions which describe the true flow in the 
whole channel (including the region near x = 0). We now extend Eq. (|2"2"j) to all x € (0,1). As explained in the 
preceding section, Eqs. ©, (fTTj) . and (f2"2"j). considered on the interval x € (0,1), define continuation of the flow 
variables to the whole physical domain. It remains only to infer an additional initial condition for the continued 
functions f'(x), U—(x), to account for the rise of differential order in the transition from Eq. (|2"Tj) to Eq. ([22)1 . To this 
end, we note that Eq. (|2ip is still satisfied by the continued functions, if it is understood that the integral term and 
the constant 7 are still evaluated on the true functions. Letting x — > then gives 2u'_(0 + ) — a/'(0 + ) = or, using 

Eq. cm, 



Wq 



a 
2 



(23) 
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This result deserves some commentary. The dimensionless quantity wq might in principle depend on the three 
independent dimensionless parameters of the theory - U, 9, and g (= gb/Uj, in the ordinary units). It is independent 
of U merely because the large-slope asymptotic expansion we use is, in fact, an expansion in powers of U, whereas 
w = 0(1). Next, independence of g is accidental: There is no reason to expect that it will hold in the subleading order 
[independence of g would be exact if wq were determined entirely by the flow properties in the narrow region near 
x = 0, since gravity is negligible on small scales (as is evident from the fact that gravitational effects are described 
by terms of the lowest differential order, Cf. e.g., Eq. (|18l) ~). However, wq is the boundary value of velocity continued 
from within the physical domain using the differential equation which involves g]. What is remarkable is that wq 
turned out to be independent of 7, and hence, of the true flow specifics near the walls. This independence makes the 
present consideration self-contained, in that it does not require the knowledge of 7 to describe flame structure far from 
the walls. This conclusion is not actually quite unexpected, as we deal here with flames evolving spontaneously, in 
contrast to flames with a fixed ignition point, where large gradients of the gas velocity near flame holder do contribute 
to wo, affecting thereby the whole flame pattern. Regarding the constant 7 itself, it should be stressed that its exact 
value cannot be found using the continued functions, precisely because they fail to describe the true flow near x = 0. 
Nevertheless, approximate values can be instructive, Cf. Sees. IIV C| tVl 
Finally, combining Eqs. (JSJ), (fTTj) one finds 

x , (1 — x)u' 

/' = u_-(l-xu'_, w-= f'\, - (24) 

u_ — (1 — x)u_ 

Substitution of these expressions into Eq. ([^1) yields a second order differential equation for U-(x), which is linear 
with respect to u'!_{x). Therefore, its solution is unique, if exists. The initial conditions for this equation follow from 
Eqs. 0, (H3, and ((231) 

u-(Q)=U, u'_(0) = ^tf. (25) 

At last, it is necessary to make the following important observation regarding general structure of equations derived 
in the leading order of the large-slope expansion. Dimensional analysis suggests that the flame speed can be written 
as 

U = C(6,gb/U})^b, 

where C is a function of the two dimensionless combinations of parameters characterizing the problem under consid- 
eration. Since the flame speed grows with g (for it is gravity that forces flame to spread along the channel), taking 
the large-slope limit implies that C(9,gb/U^) is replaced by C(9, 00), provided that the limit of C exists. That the 
limit should exist is suggested by the observation that g disappears from Eqs. (|22|) . ([24| on replacing 

However, this suffices to prove existence of the limit only if the function C (o, gb/U^j is single-valued. The nontrivial 

fact established below is that this function actually is not single- valued, and this leads to interesting complications in 
the relationship between flame speed and tube diameter, to be discussed in Sec. IIVDI 



3. Numerical solutions 



It turns out that for each set of parameters 9, g, infinitely many solutions to the equation for u-(x) exist, which 
differ by the value of U. Namely, there are two domains on the line U > 0, one of which is a finite interval U £ (Uq, Ui), 
and the other is an infinite interval U > U2 > U\, where f/0,1,2 depend on 9, g. Henceforth, solutions having U in 
these domains will be referred to as Type / and Type //solutions, respectively. Physical solutions of each type will be 
identified in Sec. IIV C[ meanwhile some examples are given in order to characterize the two types of flame patterns. 
Figure [2] shows front positions and on-shell fresh gas velocities for Type I flames with 9 = 5 and 9 = 4.5, propagating 
with a normal speed Uf = lOcm/s in a channel of b = 40 cm. The dimensionless speed of uniform movement in 
both cases is U = 11.28. Figure shows similar data for Type II flames with 9 = 7.5 and 9 = 5, both having 
Uf = 40cm/s, U = 3.9, and propagating in a channel of b = 30 cm. It is seen that in all cases longitudinal velocity 
monotonically increases with x, which implies that the gas pressure drops downstream, as it normally does in actual 
experiments. The key difference between Type I and Type II solutions is in the behavior of transversal velocity: it 
changes monotonically across the channel in the latter case, and has a pronounced maximum in the former. This 
maximum corresponds to the inflection point on the flame front, which is characteristic of Type I solutions. 
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III. FLAME PROPAGATION IN WIDE CHANNELS 

Experiments 043 show that the steady regime of flame propagation is not realized in sufficiently wide horizontal 
tubes. Instead, flames develop irregularly varying small-scale structures observed visually as a flame turbulence, or 
as a widening of the flame front, on photographs. For methane-air flames, for instance, the critical tube diameter 
d rsi 10 cm. Mason and Wheeler Q described this turbulence as "a swirling motion in a direction nearly normal to 
the direction of translation of the flame-front," and attributed it to gas convection. Figure 2] reproduces moving-film 
photographs of flames propagating in a 10% methane-air mixture in tubes of 2.5 and 20 cm diameter [6J. It shows a 
straight-line front in the former case, and complex geometrical patterns present at every instant in the latter. It is 
very important, however, that the envelope surface (the average front position) retains its shape during the phase of 
uniform movement. Moreover, it is seen from the photograph that the speed of uniform movement remains grossly 
constant. These facts justify the very notion of uniform movement, and suggest that this phase of flame evolution may 
be described in terms of flow variables averaged somehow over the small-scale flame irregularities. An appropriate 
averaging procedure can indeed be introduced, which allows extension of the obtained results to wide channels at 
the cost of introducing only one additional phenomenological parameter - the effective normal front speed describing 
amplification of the local consumption rate, caused by the small-scale front wrinkling. This will be done below 
following the lines of [19| . with several clarifications and a simplification of the averaging procedure which will be 
purely spatial (Sec. IIII A|) . The main consideration consists of two parts dealing separately with the bulk equations 
fSec. IIIIBj) . and the local propagation law (jlll C| . 



A. Averaging procedure 

Dynamics of flame disturbances characterized by the length scale A c is governed by the finite front-thickness effects. 
These effects are the result of redistribution of heat and mass fluxes in the preheat zone of the flame front, caused by 
the front curvature and inhomogeneity of the tangential flow. The simple idea behind the averaging procedure to be 
introduced is that the strong local effects due to this redistribution must largely compensate each other when their 
influence on the flame dynamics at distances 3> A c is considered. Using the condition A c <C 6, let us choose a length 
L satisfying 

A C <L<6. (26) 

Given a function A(x,t), we define its spatial average over the region V = {x : x E (x,x + Ax), As., = L} near a 
point x as 

{A) = ~hJ d2x A{x ' t] = Mx ' t} ■ (27) 

V 

By definition, (A) varies noticeably over distances ~ 6. The function A(x,t) thus turns out to be decomposed into 
two parts corresponding to the two scales b and A c : 

A = A + A 1 , (A 1 )=0. (28) 

Both the slowly varying (Aq) and rapidly varying (Ai) parts will be treated as power series in the small ratio X c /b = e 

A = 4°) + sA^ + ..., Ai = 4°> + eA^ + . . . , 

where dots denote terms of higher order in e. Despite formal resemblance of these expansions, their meaning is 
different. According to the above definition, dependence of Aq on A c is naturally expected to be regular enough 
to admit expansion in powers of X c /b. This is because A c is not the length scale of Aq, and decreasing A c only 
improves conditions (|26[) underlying the definition of Aq. Therefore, replacing A c by (eb) wherever it appears in Aq, 
and expanding in e, we obtain a series with the coefficient functions independent of A c . On the other hand, the 
principal length scale of the rapidly varying part A\ is A c , whereas the relatively weak dependence on b is expected to 
be amenable to power expansion in 1/6 (increasing b improves conditions (|26[) ). Replacing 1/6 — » e/A c therein yields 

a scries in e whose coefficient functions depend only on A c . Thus, regarding the coordinate dependence of ^4q {, it will 
be assumed in what follows that 

A r) = A r \x/b), A^ = A^(x/X c ) r = 0,l,... 
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Now, definition ([28]) implies that (A{ 0) ) + e{A^) H =0, and therefore, 

(4°)) = (4 1 )} = ...=0. (29) 

Within the approximation used throughout this paper, evolution of the flame front envelope is described in terms of 
t4q ' . This evolution is affected by the small-scale structure via nonlinear coupling of A± and in the governing 
equations. 

In subsequent calculations, the natural units are used (see Sec. Ill A[) . This implies, in particular, that A c is equal 
numerically to e, and that time t is measured in units of b/U f. We then have the following order-of-magnitude estimates 
for the functions involved with respect to the parameter e [p denotes gas pressure) 

wo = 0(1), »i = 0(l), po = 0(l), pi = 0(l), (30) 

Z,- ^ < 31 » 

d "» O(i), >g-0(i\ &-o(i\ (32) 



dxk V £ ) ®t V £ ) ' ® x k V £ 

d2V0 * 0(1), P^ = o(±), *, M = 1,2. (33) 



dxkdxi dx k dxi 
Let us proceed to the derivation of effective equations for the functions , p^ . 



B. Decoupling of dynamical equations 



The gas velocity and pressure obey the following bulk equations 

div v = , 



W + („.V)« 



-~Vp 

p 



g + ePrAv . 



(34) 
(35) 



where Pr is the Prandtl number (the ratio of viscous and thermal diffusivities, which for gases is of order unity). 
The role of the viscous term will be discussed in detail later (see Sees. UVUllV)) ; it is included already because we are 
presently concerned with the effect of dissipative processes on the large-scale flame evolution. Extracting the leading 
0(l/e)-terms in these equations with the help of estimates ([3U |) - (j3"3")) yields 



(0) 



dt 



.(0) 



,(°) 



dwvf ] = 0, 

V) v{ 0) = --Vp[ 0) 



ePrAv 



(0) 



after which collecting O(l) terms gives 



dv. 



(o) 



dt 



dt 



,(0) 



,(°) 



divw$ 0) = 0, 
V) + e ( 



,(0) 



^Vpf +g + e I Pr/\v 1 



,(o) 



,(!) 



(1) 



(36) 
(37) 

(38) 

(39) 



Equation f)38[) already involves only while the other equation couples the slowly and rapidly varying functions. 
However, the latter disappear upon averaging Eq. (l39l) . We note, first of all, that 



(i)' 



d(v\ 



(i; 



dt 



at 



= o, 



by virtue of the definitions (|2"T|) ~ (l2l?l) . Similarly, 



((^•v).r)=((^)).v) 



•V)^ o) = 0. 
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Next, 

e(VpP) = ^ 2 J dsp?, 

s 

where S is the surface of V, ds being its element. Since p^P = 0(1), the integral here is 0(L). Therefore, according 
to the choice of L, 

£ ( v ^ 1) > = (t) =o(1) - 

The same argument applies to (fg°' ) • V)^ 1 "*, as well as to the last term on the right-hand side of Eq. (l39l) . Finally, 
contribution of the last term on its left-hand side also is o(l). Indeed, integrating by parts and taking into account 
Eq. ([55]) . one has 

((„<°> ■ V) „«) = Jd*x („f • V) v?> = ^f(ds. v{V) . 
V s 

In view of Eq. (f5U|) . the last expression is 0(1/L). Hence, 

£ ((^-v)^) = o(^)=o(l). 

Thus, up to o(l)-terms, Eq. ([55)1 reduces upon averaging to the ordinary Euler equation for the functions Vq , Pq -* 

(^•v)^=-iv^+ 9 . (40) 



0t 



Thus, to the leading order, the large-scale flow dynamics in the bulk is unaffected by the small-scale wrinkling. This 
is not unexpected, of course, as the hydrodynamic equations themselves are derived using macroscopic averaging of 
the more fundamental kinetic equation. That the large-scale flow dynamics turns out to be ideal is just because we 
average precisely over the distances where the viscous forces are only important. We return to this point in Sec. lIV Cl 
The obtained result does not mean that the large-scale flame evolution is completely unaffected by the wrinkling, 
for this evolution depends also on the local propagation law of the flame front. 

C. Renormalization of local propagation law 

Account of the transport processes inside the flame front modifies the local propagation law. For the sake of 
simplicity, we shall adopt the simplest form thereof 

(v--n)-^ = l-l M ^, (41) 

where Im = 0(X C ) is the Markstein length (22J. In the gen eral case, additional terms involving time derivatives 
appear on the right-hand side of the evolution equation p3l |24| , which require modification of the auxiliary averaging 
procedure along the front, to be introduced presently. This procedure is needed to reduce Eq. (|41[) to an equation 
for the slowly varying part of gas velocity. Consider a quantity A defined at the flame front. Given a point x, choose 
Ax = Ax(x) such that the front length L between the points (x, fix, t)) and (5, f(x + Ax, t)) satisfies 

A c < £ < b, £ = 0(L). 

This is always possible since L is of the order of distance between the two points. Define the average of A over 
{x,y : x G (x,x + Ax),y = f(x,t)} by 

x+Ax 

( A )l = ^ [ dlA, (42) 
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where dl is the front line element, dl = Ndx. 

Using this operation, the unit vector n defined at the front, as well as the flame front position itself can be 
decomposed into two parts corresponding to the scales b and A c : 

n = n + n 1 , / = / + A , (m); = 0, (/i); = 0, (43) 

which can be further expanded in powers of e, as before, 

n ,i = + 0(e) , /„ = + 0(e) fx = 0(e). (44) 

Let us begin with averaging of the first term on the left-hand side of the evolution equation. In view of Eq. (|34[) . 
one can introduce the stream function ip — ?p(x, y, t) according to 

dip dtp 

U= fe' W = ~^- (45) 
Using the operation of bulk averaging (p?7|) . the function i[> can be decomposed as 

4j = ?Po + </>i , (V>i) = . (46) 

It follows from Eq. (|30j) that 

^ = O(l), ^ = 0(e). (47) 
Insertion of the decomposition (|46[) into Eq. (1431) , followed by the bulk averaging gives 



/d(^ +i>x)\ _ c^o , 1 



- i t x-\-L 



V 

and a similar equation for wq. One thus finds, up to o(l)-terms, 

Uo = ^' 100 = ~W (48) 

Next, writing n = (—f r /N, 1/-/V), and nsing Eqs. (|45|). f|47)l one has for the average normal velocity 

k+Ak x-\-Ax 



u L J \dx dzdx)_ L J 



_ y>o-|| +A - 



•o(l). 



Since tpo is a slowly varying function of spatial coordinates, the last expression can be evaluated up to o(l)-terms as 

ipo-\l + x = ipo-(x + Ax,fo(x + Ax,t),t)-il} -(x,fo(x,t),t) 

d^ d^ df \ ( df \. 

~^ h "5 — Air = u _ - w - Ax , 

cte ay ax J _ \ ox J 



where Eqs. (|44|) . ([48)1 have been used. Note that 



- df/dx\ 1 ^T^, /d/o dfA Axdf 



and similarly, 

Aa; 

n y0 = — +o(l) . 

Hence, the result of averaging the normal fresh gas velocity along the front takes the form 

(0- = K- •n )+o(l). 
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It is to be noted that no is not a unit vector. Up to o(l)-terms, its inverse norm is 



n \\- L = ; =0>1. (49) 



ox 

This quantity has a clear geometrical meaning: (/3 — 1) represents the relative increase of the flame front length due 
to its small-scale wrinkling. In general, the dimensionless quantity /3 is a slowly varying function of x. However, in the 
case of gravity-driven flame evolution we deal with (gb/Uj > 1), /3 turns out to be position-independent in zero order 
approximation with respect to e. This is because there is no appropriate length parameter to form a dimensionless 
combination involving x, on which /? might depend. Indeed, A c is not a candidate, by the very definition of f3. On 

the other hand, since /3 is a geometrical invariant, it is independent of f^°\ To prove the latter statement, we note 
that the coordinate system can always be chosen so that dfo/dx = at any given point i, in which case one can take 
Ax — L. Then dfo(x, t)/dx = O(L) for x € [x,x + L], and j3 takes the form 



-L 

dxn 



Ax L J \ y dx 



on 



(0)\ 2 



O(L), 



which proves the asserted independence of fjf^ up to terms o(l). Thus, j3 is a functional only of f[°\ and hence, is 
independent of b (Cf. Sec. MI A"|) . In other words, f3 does not vary on the length scale b and, ipso facto, on the smaller 
scale U'j/g. 

Averaging of the remaining terms in Eq. (|4ip is straightforward. Taking into account Eq. (|44[) . one finds 



x+Ax 



n) = H / dx {fo(x,t) + h(x,t)} = ^MM + o(^ 



x+Ax 



Similarly, 



X 

since /' = 0(1). Thus, the evolution equation reduces to 

(v - ■ n ) - riyofo = 1 + o(l). 
In terms of the zero-order coefficient functions, this can be finally rewritten as 



(^o° ) ^)- Z 5 r =/3: (50) 



where N = \/l + (df^/dx) , and v is the unit vector f— [/q°^'/N, ■ This result shows that from the large 



scale point of view, the effect of flame wrinkling amounts to renormalization of the normal flame speed, which in the 
ordinary units now equals to (3V /. That the renormalizing factor f3 is expressed entirely in geometrical terms, with 
no trace of the parameter Im characterizing the inner-front properties, is in accord with the idea about compensation 
of the local effects, mentioned in Sec. MI Al The estimates given in [Tj for the front area of 3D flames in wide tubes 
imply that ^ w 1.4 — 1.5. 

Finally, it follows from Eqs. (|40l) . (|50[) that jump conditions at the front for the slowly varying part of gas velocity 
are the same as in the case of zero-thickness flames. More precisely, if the velocity unit is redefined to be /3Uf, then 

the velocity jumps are given by the last two equations (|17l) upon replacing w, w, f — > m °', w °\ f^- For the sake of 
notation simplicity, one and the same designation of the flow variables will be used in channels of any width, with 
the understanding that the above replacement and redefinition of the velocity unit are made when considering flame 
propagation in wide channels. 



14 



IV. ANALYSIS OF EXPERIMENTAL DATA 
A. Relevance of 2D picture to flame propagation in tubes 

The question as to what extent the results of 2D analysis are applicable in 3D situation is quite nontrivial, and 
in each particular problem requires special consideration. To settle it in the present case, we turn to Fig. [5] which is 
a series of snap-shot photographs of the early stage of flame evolution after ignition, made at equal time intervals, 
in a glass tube of d = 10 cm |6j. It shows that as the flame travels along the tube, its initially hemispherical front 
not only tilts but also flattens, as is evident from the fact that the front image becomes thinner. This observation 
suggests that the flame shape in the regime of uniform movement can be described using 2D tools, and implies that 
the channel width is to be set equal to the tube diameter. As to the speed of uniform movement, it was proposed in 
[ill l25l | that the 2D values are to be multiplied by a factor w 1.5 to obtain the flame speed in a tube (of diameter 
equal to the channel width), the reason being that this relation holds for bubbles of similar shape. However, as was 
already mentioned in Introduction, the theory of bubble motion is irrelevant to flame propagation. As a matter of 
fact, there is a solid reason to expect that no such factor is to be introduced in flame theory. The point is that unlike 
bubbles, flames propagate because they consume fuel; the rate of this consumption is proportional to the flame front 
area, which is determined by the slope of front elements, rather than their azimuth with respect to the tube axis. 
Thus, the relative velocity increase of 2D and 3D flames is the same if they have the same average slope. 

The misleading analogy with bubbles helps to consider this issue from a somewhat different standpoint. To get 
insight into how flame dimensionality might affect its speed, compare propagation of flames and bubbles in a channel 
and in a tube of equal size. To simplify geometry, it is more convenient to deal with vertical propagation, wherein 3D 
shapes can be assumed axisymmctric. Consider first bubble motion, and let the 3D bubble shape be obtained from 
the corresponding 2D channel pattern by rotating it around the tube axis. The bubble speed depends essentially on 
the relative cross-sectional area of the jet formed between the bubble and the walls, which is evidently larger in the 
3D case. The same geometrical reason suggests that if it were possible to obtain an axisymmetric flame by rotating 
a 2D pattern around the tube axis, then the speed of this flame would indeed be larger than the speed of its 2D 
original as proposed in [TP. l25j]. since the front slope grows with distance from the tube axis. However, such rotation 
would violate flow continuity, for the gas elements burnt near the tube wall come from regions of lesser radial distance 
from the axis (Cf. Fig. [SJ, whereas the flow far upstream remains homogeneous. Therefore, instead of increasing 
the speed of propagation, 3D flame has to reduce its longitudinal spread to maintain flow continuity. On the other 
hand, retaining continuity in the case of bubbles of fixed shape is not a problem, for their motion is not bound to 
satisfy such condition as the constancy of local consumption rate. We thus see that the difference in the mechanisms 
governing propagation of bubbles and flames implies essentially different sensitivity to the flow dimensionality. 

It can be added that in complete analogy with the 2D case, the effect of small-scale 3D front structures amounts to 
renormalization of the local propagation law (the proof is almost identical to that given in the preceding section). To 
be sure, these arguments do not affirm equality of propagation speeds of 2D and 3D flames, but they point toward 
certain rigidity of the flame speed with respect to manifestations of the three-dimensional nature of the problem, such 
as occurrence of pronounced 3D front structures, or circularity of the tube cross-section. On these grounds, the speed 
of uniform movement of flame in a tube of diameter d will be identified henceforth with the speed of steady 2D flame 
propagating in a channel of equal width, b = d. 

B. Assessment of heat losses 

For flat-on-average flames, heat losses to the tube walls are usually of minor importance in regard to the flame speed, 
for their relative value rapidly decreases with increasing tube diameter. Thus, although preventing flame propagation 
in very narrow tubes, heat losses are practically negligible in tubes of a few centimeters diameter. For instance, the 
speed of downward propagation of planar methane-air flames near the limits of inflammability (in which case the 
gravity and transport effects completely stabilize planar front) is the same in tubes of d = 5 cm and d — 23 cm p6| . 
implying that even for slowest flames the cooling effect is negligible already in tubes of 5 cm diameter. Things are 
different, however, for horizontally propagating flames, because the large flame spread along the tube considerably 
enhances heat outflow from the front, reducing thereby the flame speed. The rate of this outflow is, in the ordinary 
units, 

q = Jf EAT, 

where E s» ird x (U / j3U f)d/2 is the area of the tube surface surrounding the front (the factor 1/2 accounts for the fact 
that heat is transferred only from the hot side of the front), AT the temperature difference between the cold wall and 
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burnt gas, and K the heat transfer coefficient. On the other hand, the longitudinal heat flux due to deflagration is 

q = pQ xU x 7rd 2 /4, 

where Q is the specific heat of combustion, and p is the cold gas density. The dimensionless ratio 

q 2K 

8 =- q = i^m> (51) 

where c p = Q/AT is the average specific heat capacity of fresh gas, gives the relative velocity drop due to the heat 
losses to the walls. It is seen that this parameter depends neither on d, nor on U. In application to flame propagation 
in methane-air mixtures under standard conditions, p 1.2kg/m 3 , c p ranges 1.3 to 1.8 kJ/kg-K, depending on 
the methane concentration, while the heat transfer coefficient for tubular air/metal heat exchangers ranges 5 to 
35W/m 2 -K [27l ]. As no relevant information is given in the cited works 5 to further specify this parameter, the mean 
value K = 20 W/m 2 -K will be used in all subsequent estimations of 8. Substitution of the above figures into Eq. ([ST]) 
shows that the velocity drop can be as large as 30% (this value is attained near the limits of inflammability) . 



C. Identification of physical solutions 

The family of all possible steady flame patterns was described in Sec. Ill B 31 To select practically relevant solutions, 
one has to invoke properties of the physical gas flow existing far downstream. The point is that despite ideality of 
the large-scale flow, proved in Sec. MI B[ its asymptotic structure downstream is essentially determined by the gas 
viscosity. Indeed, even if viscosity is so small that its effect is negligible at finite distances from the front, it leads 
to gradual decay of vorticity produced in the flame front, whereby the flow becomes homogeneous sufficiently far 
downstream (for simplicity, we discard viscous drag exerted by the walls, which results in the logarithmic velocity 
profile) . It is worth recalling in this connection that the master equation ((4} itself was derived in the sense of analytic 
continuation to the limit of zero viscosity fl8j . 

To identify relevant solutions, we note that the vorticity decay is accompanied by a loss of kinetic energy of the 
burnt gases. Hence, criterion for selecting physical solutions follows from the energy conservation law: the total loss 
of kinetic energy downstream must be compensated by its gain in the process of combustion, i.e., by the gravitational 
work done on gases in the transition domain. 6 The latter is the region where transversal component of the gas velocity 
is nonzero (the gravitational power vanishes together with w). The reasoning used to derive Eq. (|14p implies that 
this is the region indented in Fig. [5] by two vertical broken lines. The above criterion can be expressed in terms of 
the kinetic energy flux through a vertical section at the point y, to be denoted qr{y)'- the total kinetic energy loss 
downstream is the difference qrioo) — qriU), whereas the kinetic energy gain is qr(U) — (?t(0). Equating the two 
expressions yields 

qT(U) = ~{q T (oo)+q T (0)} . (52) 

By virtue of the total mass conservation, velocity of the asymptotic homogeneous flow is, in the natural units, 8U 
(recall that according to our definition of the velocity unit in wide channels, this is 9U[3Uf, in the ordinary units). 
Hence, 

^{q T (oc) + q T (0)} = - \eU y —f + U— \= K - -I—. (53) 



2 2 ^ 2(9 2 

To express the left-hand side of Eq. (|52j) in terms of the fresh gas velocity, we note that the flow continuity and 
Bernoulli integral imply the following relations between the on-shell burnt gas velocity and pressure, and their values 



5 That the heat losses to the tube walls can significantly affect flame propagation was recognized already by Mallard and Le Chatelier 
who estimated them for flat flames and concluded that the effect is negligible in sufficiently wide tubes, the critical tube diameter being 
approximately 5 cm (Cf. Q). This conclusion was apparently shared by Coward and Hartwell despite their discovery that the large 
increase of flame speed in horizontal tubes is the result of large longitudinal flame spread. 

6 In the leading order of the large-slope expansion, the kinetic energy increase due to velocity jump across the front is negligible in 
comparison with the gravitational work. 
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at the line y — U: 

0N(rj)drj = u{x,U)dx, 

. . . u 2 (x,U) 
-J- h 6p + {r)) - gr\ = h 6p{x, U) - gx, 

where (x, U) is the intersection of the line y — U with the stream line that crosses the front at (77, f(r))), Fig. [5J Since 
w(x, U) — 0, the combination [9p(x, U) — gx] is independent of x, so that it can be replaced by [9p+(l) — g\. Taking 
into account also that p+(l) = P-(l) — a, and using Bernoulli integral for the upstream flow, we find 



q T (U) 



dxu(x, U) 



u 2 (x,U) 
20 



1 

- J drjf'irj) 



u 2 _{l)-U 2 



20 



-g(v- !) 



With the help of the relations /(l) — U, f = U— — (1 — x)u'_, this can be finally written as 



<ZT(t/)^K(l)-C/ 2 ]+ V 



drj [% 



(1 - T,)U'_] 



(54) 



The right-hand sides of Eqs. (|53l) . (|54p are functions of U. Physical solutions correspond to intersections of their 
graphs, as illustrated by Fig. [7] in the case 9 — 7, g — 61.3. The left solid curve is the plot of qr(U) for Type I 
solutions U € (4.45,4.753), whereas the right solid curve represents the initial part of qr(U) for Type II solutions 
U > 6.615; the dotted curve is the plot of (9 2 + l)Z7 3 /4. It is seen that there is one physical Type I solution with 
U = 4.72, and one physical Type II solution with U = 6.96. Numerical scrutiny of the equation for u_ has shown that 
this situation is generic: For every set of parameters 9, g, there is exactly one physical solution of each type. 

To conclude this section, let us return to the issue touched upon at the end of Sec. Ill Al namely anticipated vanishing 
of the last term in Eq. Q. As was already noted in Sec. Ill B 21 its exact value cannot be found using the functions we 
work with, because they are obtained as solutions to Eqs. JSJ, (fTTj) . ([22]) continued over the narrow region near x = 
where the large-slope approximation is not valid. Despite its narrowness, this region gives rise to a non-negligible 
contribution to the integral (|20j) . for the integrand is proportional to the derivative of gas velocity which rapidly varies 
near the wall. Under such circumstances, making use of approximate values cannot decide whether 7 vanishes on 
the physical solutions, but it gives valuable qualitative information about its behavior on the family of all solutions. 
Specifically, combining Eqs. (IT7|) . (|T8)) . the constant 7 can be written as 



drjw^ 



din v + 
drj 



9_ 



Taking into account that w+ ~ a, this expression shows that in general, 7 ~ a 2 . For instance, in the above example 
7 = 52.4 for U — 4.6 (Type I solution), and 7 = 32.1 for U — 8.0 (Type II solution). However, the approximate value 
of 7 significantly drops down on the physical solutions: its values for Type I and Type II physical solutions are 9.2 
and 13.5, respectively. Moreover, 7 has a root in a vicinity of the physical solution: 7 = for U = 4.735 and U — 6.76. 
This observation gives support to the hypothesis that the exact 7 vanishes on physical solutions. If true, this would 
deliver an intrinsic criterion for selecting physical solutions, without the need to invoke asymptotic properties of the 
flow downstream. 



1. The range of applicability of the theory. Accuracy assessment 



The asymptotic formula ([S]) for M-operator requires for its validity only that the front slope be large, /' 3> 1, or 
equivalently, [/> 1, and so evidently does Eq. ((HJ. The relative value of neglected terms in these formulas is 1/U 2 . 
However, derivation of Eqs. (jlip . (|2"21 essentially uses the estimates u = 0(U), w — 0(1) to omit w in comparison 
with u. Taking into account that the magnitude of w is a (since wq = a/2), this implies a stronger condition U 3> a. 
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Switching to the ordinary units and recalling that U ~ \fgd , we thus obtain the formal criterion of applicability of 
the theory 

aUf 

1 < 1. 



vgd 

The meaning of this condition is that the ratio on its left-hand side estimates the relative error introduced by the 
above-mentioned omission. 

Even without detailed analysis of solutions, it is clear that flame propagation with U < a/2 is not covered by the 
present theory, since Eqs. (fT4")l . (|2"3")l imply that in this case u_(0) < tu_(0), and so the relative error exceeds unity, 
at least at the flame head. Hence, solutions obtained in Sec. Ill Bl accurately describe the front position and velocity 
profiles, only if U a/2. As to determination of the physical value of U itself, the requirement is not actually that 
stringent, for U is an integral quantity. The dominant contribution to U comes from the region where the front 
slope is largest, that is, from the lower part of the channel. In Type II solutions, U-{x) increases with x, u'_ being 
proportional to a, while W-(x) monotonically decreases, so that condition u_ ^> w_ can be considered practically 
satisfied, provided that it_ > w_ at the upper wall, i.e. if 

U>~. (55) 

In Type I solutions, the u-component of gas velocity behaves as before, while the transversal component first increases, 
picks at x ~ 0.2, where it is generally of order w_, and then rapidly decreases. Since the pick is in the upper part of 
the channel, criterion (|55p can be applied to this type of solutions as well. Accuracy of the theoretical value U within 
the range can now be estimated as follows. Let r denote the relative error in U. This error comes mainly from the 
region where w/u > r. On the other hand, since U is equal to the front length, r is approximately the length of the 
front arc belonging to this region, relative to the total front length. Therefore, r can be determined as the ordinate 
of intersection of the curves w/u and f /U. An example is given in Fig. [8] for Type I solution with 6 — 5, b — 40 cm, 
Uf = lOcm/s, and U = 11.3, in which case r 10%. The relative error is smaller in wider channels and for slower 
flames, being about 20% at the boundary of the domain defined by ([55]). As subsequent analysis shows, larger errors 
outside this domain make comparison of the theory and experiment inconclusive. 



D. Comparison with experimental data on methane-air flames 



We now go over to a detailed comparison of the obtained results with the experimental data on methane-air flames, 
given in . The flame speed data to be analyzed is collected in figures 2 and 3 of [gj], the first of which is reproduced 
in Fig. [H] The experimental marks are due to Mason and Wheeler (open circles) Q, Mason (triangles) Q, Coward and 
Greenwald (crosses) [H, and Coward and Hartwell (filled circles) Q. The curves represent an attempt to interpolate 
the experimental data for various methane concentrations. Inflections in the curves in the region 10 cm< d < 20 cm 
indicate the onset of formation of complex front structures observed in wide tubes. 

To begin with, we note a very important evidence contained in the table reproduced in Fig. [9j namely that the 
experimental data for the speed of uniform movement (denoted S.U.M. in the table) cannot be approximated by a 
single function of the form 

U = C{6)^g~d, (56) 

suggested by the naive dimensional analysis. To be more specific, consider a pair of flames in a tube of 20 cm diameter, 
corresponding to c = 8.02% and c = 12.05%. According to the table, they propagate with the speeds U — 102cm/s 
and U = 84 cm/s, respectively. On the other hand, both flames have ss 6.9 [3(|. Furthermore, the heat loss 
parameter 6 also has the same value 0.07 in the two instances. Therefore, according to the naive formula, the flames 
ought to travel with the same speed. A similar pair for d = 10 cm is c = 8.21% and c = 11.94%, in which case 
6 = 6.94, and 5 — 0.064, but the corresponding flame speeds are 84 cm/s and 68 cm/s. That the differences in the 
observed values are well above experimental errors is clearly seen from the same table which contains data for pairs of 
close concentrations at the points being compared. For instance, the above figure for c = 8.02% comes together with 
that for c = 8.15% (U = 103cm/s), whereas the one for c = 12.05% - with c = 12.04% (U = 81 cm/s). Apparently, 
the authors of 6] tabulated the data so carefully because they recognized the difficulty of fitting the flame speed data 
into ansatz (|56|) . though did not state this explicitly. 7 



7 These authors emphasize that the relation between flame speed and tube diameter cannot be expressed in the form U = Cd k with 
constant C, k, but relate this to the existence of inflections in the curves U(d). 
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Directly related to this is a peculiarity exhibited by the flame speed considered as a function of methane concentra- 
tion, U(c), namely, existence of inflection points at c ~ 8% and c pa 11.4%. It is seen from the data in Fig.[3]that after 
a relatively moderate rise between c = 6% and c pa 7.2%, U(c) increases significantly more rapidly on the interval 
c G (7.2%, 9%), after which it remains almost constant for c up to 10.4%. Similarly, the fall of U(c) on the interval 
c G (11%, 12.2%) is considerably larger than on (12.2%, 13.5%). The ratio of the flame speed at the central plateau 
to its values at the foothills (i.e., at c pa 7.2% and c pa 12.2%) exceeds 1.5. This behavior cannot be explained by 
differences in the normal flame speed at different fuel concentrations, for U is independent of Uf when gd > Uf. 

Proceeding to explanation of the observed flame behavior, let us first identify the factors that determine variation 
of the flame speed with fuel concentration. They include 

(1) Dependence of the flame speed on the gas expansion coefficient, 

(2) Heat losses to the tube walls, 

(3) Existence of two distinct regimes of steady flame propagation. 

Regarding (1), numerical analysis of the differential equation for U—(x) shows that the flame speed of physical Type I 
(Type II) solutions increases (decreases) as 9 increases, though this dependence is relatively weak. Next, heat losses 
significantly affect propagation of flames near the inflammability limits, giving rise only to a few-percent correction 
to the flame speed in near-stoichiometric mixtures. At last, the existence of two types of solutions with the flame 
speed ratio about 1.5 is the key factor for explaining the aforementioned peculiarities in the flame speed data. As 
will be shown below, the inflections in U(c) correspond to transitions between Type I and Type II solutions. The 
qualitative picture of the dependence U(c) is thus as follows. Mixtures near the limits of inflammability and near- 
stoichiometric mixtures sustain respectively Type I and Type II regimes of flame propagation. The joint effect of 
factors (1),(2) explains the moderate rise of the flame speed on the interval c G (6%, 7.2%), and also its moderate fall 
on (12.2%, 14%). In near-stoichiometric mixtures, on the contrary, these factors counteract partially compensating 
each other, which produces the plateau in U(c) for c G (9%, 10.4%). Transitions between the two regimes explain the 
rapid changes of the flame speed in mixtures with c G (7.2%, 9%) and c G (11%, 12.2%). 

To carry out quantitative comparison of the theory and experiment requires some auxiliary experimental data. The 
normal velocities deduced in Q from the front area calculations are somewhat lower than the planar flame speeds, 
because of discard of the finite front-thickness effects. The actual values of Uf will be adopted from [28, 29]. The gas 
expansion coefficients are based on the adiabatic flame temperature calculated with the help of [3(|. The heat losses 
are taken into account by multiplying the theoretical value U by (1 — 8), where 6 is given by Eq. (|51[) . In this formula, 
p is set equal to 1.2kg/m 3 (neglecting its variation with methane concentration) and K = 20W/m 2 -K in all cases 
considered below. The quantity c p is calculated in each instance using the data provided in [3(|. Regarding /3, the 
estimates given in Q imply that /§ pa 1.45 in a tube of 24 cm diameter. In the absence of more detailed information, 
f3 = 1.45 will be used in all wide tubes (d 20 cm). In tubes of moderate width, (3 — 1, by definition. 

Let us start from the lower inflammability limit represented in Fig. [5] by the lowest curve interpolating experimental 
data for c = 6%. The following is a detailed example of the computational procedure used below to compare the theory 
and experiment at various methane concentrations. Recalling that the stoichiometric methane/air composition is at 
c = 9.5%, one first finds the equivalence ratio for the mixture under consideration: <fr — 0.61. It is then used to 
determine adiabatic flame temperature with the help of [30], which for standard conditions (latm, 300 K) yields 
1680 K. Hence, AT = 1380 K and 9 — 5.64. Type I and Type II physical solutions found numerically for this 9 read 

Ui = 18.2Vd, Un = 29.3Vd, 

where d is to be expressed in centimeters. Next, to take into account the heat losses, one uses 0.9MJ/mol for 
methane's heat of combustion, together with the found AT, to compute c p pa 1.3kJ/kg-K. The normal flame speed 
Uf = llcm/s is inferred from Fig. 2 of [2{|- Substitution into Eq. (l5"Tj) gives S = 0.16 and S = 0.23 for tube 
diameters d ^ 20 cm and d ^ 10 cm, respectively. Finally, condition (|55|) implies that the above numerical solutions 
are valid as far as U > 2.3, or in the ordinary units, U > 2.3 x llcm/s pa 25cm/s in tubes of moderate width, and 
U > 2.3 x 1.45 x llcm/s ~ 40cm/s in wide tubes. It is seen from Fig.|9]that for the given methane concentration, the 
range of applicability of the theory covers all available experimental data. The computed flame speeds for Type I and 
Type II solutions, combined with this data, are plotted in Fig. ITUT a). Clearly, flames propagating in a 6% mixture 
favor regime described by Type I solutions (solid curve) in tubes of any diameter. The same is true of 7% mixture, 
Fig. [TOTb). where Uf = 20cm/s, 9 = 6.26. In this case, condition (f55|) requires U > 53cm/s, hence only results for 
d ^ 10 cm are shown. On the contrary, burning in mixtures close to the stoichiometric composition is of Type II, as 
is evident from Fig. [TT1 representing U(d) for 9% and 10% mixtures: within the error of calculation, all experimental 
marks in these two cases belong to the curve representing Type II solutions (dotted line). The computed values of 
the main parameters used to draw the curves for various mixtures are collected in Table I, together with the relevant 
parameters of physical solutions. The latter include the constants U(d)/y/d for each type of solutions, and the error 
of calculation determined as described in Sec. IIV C II Given in the last two columns are pairs of reference values of 
the error - r m for d — 10 cm, and r w for d = 100 cm. The figures for the heat loss parameter <5 refer to wide tubes; 
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in tubes of d ^ 10 cm they are to be multiplied by 1.45. Regarding mixtures with c = 8% to 12%, condition (|55|) 
implies that the theoretical results are amenable to verification only in wide tubes. Taking c = 10% for instance, one 
has Uf = 43cm/s, 9 — 7.5, and Ui = 1.4, Uu = 2.03 in a tube of 10 cm diameter, whereas a/2 — 3.25; to put it in 
terms of accuracy, the error of calculation r = 23% makes comparison of the theory and experiment inconclusive in 
this case. 

The above consideration leads us to the conclusion that the rapid change in flame speed, observed in the range of 
concentrations c = 7.2% to c = 9%, represents transition between Type I and Type II regimes of flame propagation. 
To get insight into the nature of this transition, let us turn to Fig. [T2] containing the flame speed data in 8% mixture. 
The following consideration shows that this data points toward existence of an intermittency of the two regimes in 
transient mixtures. Namely, it is seen that while the experimental marks for d = 30.5 cm are clearly of Type II, the 
mark for d = 96.5 cm belongs to neither Type I nor Type II solution. In fact, the error of calculation for this tube 
diameter is r — 12% for Type I and r = 10% for Type II solution, so the distance of the mark from cither curve is 
well above the error of calculation (this mark is due to [H according to which the experimental error is about 2-3%). 8 
This observation can be easily explained if we admit that the flame undergoes transitions between the two regimes 
as it travels in the tube. This explanation is suggested by the very fact that the given mixture is transient, and is 
most naturally formulated in terms of stability of the solutions. Indeed, that the mixtures with c < 7% propagate 
only Type I flames means evidently that Type II solutions are unstable in these mixtures (or, at least, that they are 
considerably less stable than those of Type I), and that situation is opposite in mixtures with c > 9%. It is therefore 
natural to expect that in the 8% mixture, the two types of solutions are approximately even regarding their stability 
properties. Intermittency occurs if both solutions are unstable or neutrally stable. 

The coexistence of two solutions with similar stability properties is not sufficient to establish intermittency in the 
case when they are stable, for it would not exist without a mechanism that continuously provokes switching of the 
propagation regimes. However, such mechanism is readily identified. Namely, it is not difficult to see that fluctuations 
of the heat flux due to heat losses into the walls are sufficient to trigger transition from one regime of flame propagation 
to the other. For this to happen, the heat inflow/outflow, AQ, necessary to accelerate/decelerate the flame must be 
comparable to the difference of kinetic energies of gases in the two regimes, and occur on a time interval less than the 
characteristic time of perturbation damping. For the gravity-controlled flame evolution this time is of order Uf/g. 
Multiplying it by q from Sec. HVBl yields 4 x 10 3 W for Type I solution (U = 4.25), or in the natural units, 5 x 10 4 ; 
for Type II solution (U — 6.3) the figure is 1.5 times larger. On the other hand, the kinetic energy of gases in the 
volume surrounded by £ can be evaluated as 

1 ' 29 4 ' 

which gives Ek{Ui) ~ 2 x 10 2 and Ek{Uu) ~ 7 x 10 2 . Therefore, even a few-percent local fluctuation of the heat 
transfer coefficient will produce AQ sufficient to switch the propagation regime. The resulting intermittency can in 
principle be observed as a recurring appearance of an inflection point at the front head (Cf. Figs. (2J [3]). The above 
stability consideration implies that flames in transient mixtures are particularly susceptible to such fluctuations, but 
the described mechanism may work also in mixtures where the two regimes of flame propagation are of markedly 
different stability. In that case, a more intensive fluctuation of the heat flux is required to perturb the flame, as the 
damping time is shorter in more stable regimes, and no proper intermittency occurs. 

For an experimental evidence of the proposed scenario, let us turn to the work Q remarkable, in particular, by its 
scrupulous exposition of the experimental procedures. Regarding flame propagation in a tube 30 m long and 30.5 cm 
in diameter, the authors mention that "some experiments demonstrated that there was a varying of the flame speed 
throughout the tube, an occurrence which was more frequent the more disturbed the external atmosphere," and 
emphasize that experiments were made only when the wind was negligible near the open end of the tube (Cf. pp. 17, 
20, and 22 in The latter precaution excludes direct atmospheric impact on the burnt gas flow, implying that the 
observed variations in the flame speed are caused by variations of the external conditions along the tube. That the 
value of the heat transfer coefficient strongly depends on the ambient air conditions is well known. 

Returning to consideration of the function U (c) , one finds that the flame speed behavior at methane concentrations 
larger than 10.4% is quite similar to that at c < 10.4%. Namely, comparison of the theoretical and experimental 
data for a 11% mixture, Fig. I13f a). shows that this mixture sustains Type II regime, whereas flames propagating 
in 13% and 14% mixtures are of Type I in tubes of any diameter, as is evident from Figs. ir^Tcl fTST d'). and as was 
contemplated in the qualitative analysis of the experimental data, given earlier. Figure rT3Tb) suggests that the 12% 



As to the mark for d = 20 cm, the data is inconclusive, since the two curves are close to each other, whereas calculational error is 
relatively large for this tube diameter, r = 20%. 
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mixture is transient, again in conformity with the observation that the flame speed decreases most rapidly on the 
interval c € (11%, 12.2%), although proximity to the right end-point of this interval does not give good grounds to 
affirm intermittency in this mixture: the distance of the experimental marks for d = 30.5 cm and d = 96.5 cm from 
the curve representing Type I solution only slightly exceeds the calculational errors (r = 16% and 12%, respectively). 

V. DISCUSSION AND CONCLUSIONS 

Analysis of the horizontal flame propagation, carried out in Sec. HH revealed existence of two different types of 
steady flame patterns formed under the strong influence of the gravitational field. Characteristic of Type I flames is 
the existence of an inflection point at the flame front head, with an associated pronounced maximum of transversal 
component of the on-shell fresh gas velocity. The longitudinal component of the gas velocity increases monotonically 
along the front in both cases, but more rapidly in Type II flames which always propagate faster than Type I flames 
under the same conditions. Next, it was found in Sec. IIV Cl that there is exactly one solution of each type, satisfying 
the physical requirement of asymptotic homogeneity of the burnt gas flow. At last, detailed comparison with the 
experimental data on methane-air flames, performed in Sec. IIVDI showed that both regimes are realized in practice: 
mixtures near the limits of inflammability propagate only Type I flames, whereas flames in mixtures close to the 
stoichiometric composition follow the faster Type II regime. This result resolves the difficulty mentioned in Sec. IIVD1 
namely that the experimental data for the flame speed versus tube diameter cannot be fitted into a single-valued 
function: Two flames with equal 9 have different velocities in tubes of the same diameter because they propagate in 
regimes of different type. Furthermore, inflection points in the functions U(c) correspond to transitions between the 
two regimes of flame propagation, which take place in the regions characterized by rapid changes of the flame speed. 

The comparison also clarifies the origin of infections in the curves representing flame speed versus tube diam- 
eter. Namely, the small-scale wrinkling that develops in wide tubes effectively increases the normal flame speed, 
reducing thereby the heat losses: according to Eq. (ISTl) . they are damped by a factor j3 = 1.45. This effect is espe- 
cially pronounced in mixtures near inflammability limits where the heat losses are significant, and is clearly seen in 
Figs. [T3Tc) JTST d) . Of course, this leaves open the question why complex front structures develop in tubes of d > 10 cm, 
but this is entirely a stability issue. Incidentally, occurrence of these structures cannot be described as a result of de- 
velopment of the Darrieus-Landau instability: Behavior of flame perturbations in the highly nonlinear flows generated 
by flames propagating with U 3> Uf is completely different from that found for planar flames [17| . In this connection, 
a general fact to be emphasized is that formation of steady flame patterns as well as dynamics of disturbances of 
flames with U ^> Ut is essentially determined by vorticity production in the flame. In the considered case of horizontal 
flame propagation, this process is controlled by the baroclinic effect. 

A technical issue that deserves special discussion is the identification of physical solutions. As was already mentioned 
in Sec. IIV CI the requirement of asymptotic homogeneity of the burnt gas flow, used to select physical solutions, is 
intimately related to the fact that within the on-shell approach, flow dynamics is considered in the limit of vanishing 
viscosity. It was also pointed out that the criterion identifying physical solutions may probably be formulated as the 
requirement of vanishing of the parameter 7 that describes effective gravity in the system, although this cannot be 
asserted based on approximate evaluation of 7 in Sec. IIV CI What can be affirmed, however, is that 7 is certainly 
nonzero far from physical solutions, and in view of equivalence of gravity and translatory acceleration, this has an 
important implication regarding the phenomenon of flame self-acceleration. To explain the point, let us allow for small 
but nonzero viscosity. This introduces no appreciable changes in the structure of Type I or Type II solutions, but 
significantly changes criterion for selecting physical solutions. The reason is that the viscous drag exerted by the walls, 
which now cannot be neglected, causes additional energy losses by the gas flow, increasing with the distance traveled 
by the flame. For sufficiently small viscosity (such that the viscous length > bU) this process can be considered 
quasi-steady. To compensate for the additional losses, a larger kinetic energy flux qr(U) is required. An immediate 
consequence is that this causes deceleration of Type I flames, and acceleration of Type II flames, as is evident from 
Fig. [7] Since Type II solutions exist for all U > U2, flame acceleration continues until compressibility effects come 
into play. This observation opens a way to describe the flame pre-acceleration which is well known to be an essential 
(and least studied theoretically) stage of deflagration-to-detonation transitions. 

It is to be noted, finally, that the methods developed in this paper can be used almost directly to study high-speed 
regimes of vertical flame propagation. 
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TABLE I: Parameters used to draw theoretical curves in Figs. 10-13. 
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Figure captions 



Fig.l: Flame propagation in horizontal channel of width b. (a) Symmetric pattern obtained by reflection of the 
physical channel with respect to y-axis. Vertical arrows depict acceleration of gravity, g, which in the image domain 
turns out to be directed upward, leading to a discontinuity in g at x = 0. (b) Physical realization of the discontinuity 
by a gravitating layer placed at x — 24 

Fig. 2: Front position and on-shell fresh gas velocity of Type I flames with 6 — 5 (solid line) and 9 — 4.5 (dotted 
line), propagating in a channel of b = 40cm; normal flame speed Uf = lOcm/s, dimcnsionlcss speed of uniform 
movement Z7 = 11.28 25 

Fig. 3: Same for Type II flames with 9 = 7.5 (solid line) and 9 = 5 (dotted line); b = 30cm, Uf = 40cm/s, U = 3.9. 

26 

Fig. 4: Moving-film photographs of flames propagating in 10% methane-air mixture in tubes of 2.5 cm (left) and 
20 cm (right) diameter. Both flames travel 26 cm distance. Source: Reproduced from Plate I of [6j 27 

Fig. 5: Series of snap-shot photographs of the early stage of flame evolution after ignition, made at equal time 
intervals in a glass tube of d — 10 cm. Source: Reproduced from Figure 10 of [|| 27 

Fig. 6: Schematics of steady flame propagation in horizontal channel. Horizontal arrows of different length qualita- 
tively depict gas velocity profiles in various parts of the channel. Also shown are a pair of neighboring streamlines, 
and transition domain indented by two vertical broken lines - the region where uj^O 27 

Fig. 7: Identification of physical solutions in the case = 7, g = 61.3 (natural units). Solid curves represent qr(U) 
for Type I (left) and Type II (right) solutions, dotted curve is the plot of (9 2 + l)C/ 3 /4 28 

Fig. 8: Determination of the relative error in U for Type I solution with 9 = 5, b — 40cm, Uf = lOcm/s, and 
U = 11.3. Solid curve is the plot of w/u, dotted curve - f/U. The ordinate of their intersection s» 0.11 28 

Fig. 9: Collection of experimental data for the speed of uniform movement (S.U.M) in tubes of various diameters. 
Source: Reproduced from page 1999 of [|| 29 

Fig. 10: Comparison of Type I (solid line) and Type II (dotted line) numerical solutions with the experimental data 
for flame speed in 6% (a) and 7% (b) methane-air mixtures. Theoretical curves are corrected by the factor (1 — S) 
accounting for the heat losses. Flame speed is measured in cm/s, tube diameter - in cm 30 



Fig. 11 
Fig. 12 
Fig. 13 



Same for 9% (a) and 10% (b) methane-air mixtures 30 

Same for 8% methane-air mixture 31 

Same for 11% (a), 12% (b), 13% (c), and 14% (d) methane-air mixtures 32 
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